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ABSTRACT 

This is the third paper of a series describing the Asiago-ESO/RASS QSO 
survey (AERQS), a project aimed at the construction of an all-sky statistically 
well-defined sample of relatively bright QSOs (B < 15) at z < 0.3. We present 
here the clustering analysis of the full spectroscopically identified database (392 
AGN). The clustering signal at 0.02 < z < 0.22 is detected at a 3 — 4a level and 
its amplitude is measured to be r = 8.6 ± 2.0/i -1 Mpc (in a AC DM model). 
The comparison with other classes of objects shows that low-redshift QSOs are 
clustered in a similar way to Radio Galaxies, EROs and early-type galaxies in 
general, although with a marginally smaller amplitude. The comparison with 
recent results from the 2QZ shows that the correlation function of QSOs is con- 
stant in redshift or marginally increasing toward low redshift. We discuss this 
behavior with physically motivated models, deriving interesting constraints on 
the typical mass of the dark matter halos hosting QSOs, Mdmh ~ 10 12 ' 7 /i _1 M Q 
(10 12 - - 10 13 ' 5 /i _1 M o at la confidence level). Finally, we use the clustering data 
to infer the physical properties of local AGN, obtaining Mbh ~ 2 • lO^^M© 
(10 7 - 3 • lO^^M©) for the mass of the active black holes, t agn ~ 8 • 10 6 yr 
(2 ■ 10 6 — 5 • 10 7 yr) for their life-time and r] ~ 0.14 for their efficiency (always for 
a AC DM model). 
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Subject headings: Surveys - Quasars: general - Clustering: quasar - Cosmology: 
observations 



1. Introduction 

The analysis of the statistical properties (luminosity function and clustering) of the 
cosmic structures is a fundamental cosmological tool to understand their formation and 
evolution. The clustering of QSOs and galaxies at small to intermediate scales (1-50 h^ 1 
Mpc) provides detailed information on the distribution of Dark Matter Halos (DMHs) that 
are generally thought to constitute the "tissue" on which cosmic structures form. This 
means to investigate -indirectly- fundamental astrophysical problems, such as the nature of 
dark matter, the growth of structures via gravitational instability, the primordial spectrum 
of density fluctuations and its transfer function. The lighting up of galaxies and other 
luminous objects, such as QSOs, involves complex and non-linear physics. It depends on 
how the baryons cool within the DMHs and form stars or start accreting onto the central 
black hole (BH), ending up as the only directly visible peak of a much larger, invisible 
structure. The so-called bias factor, b(r,z), is used to explain the difference between visible 
structures and invisible matter, whose gravity governs the overall evolution of clustering. 
This complex relation is summarized by the simple formula £(r, z) = b 2 (r, z)£ m (r, z), where 
£(r, z) and £ m (r,z) are the two-point correlation functions (TPCF) of radiating objects and 
dark matter, respectively. In this way the detailed analysis of the distribution of the peaks of 
visible matter can distinguish among the various models for the formation of structures. In 
particular the hierarchical growth of structures is naturally predicted in a cold dark matter 
(CDM) scenario, where larger objects are constantly formed from the assembly of smaller 
ones. An alternative view of the structure formation and evolution, supported by both 
some observed properties of high-redshift ellipticals and EROs (Daddi et al. 2001, 2002) 
and theoretical modeling (Lynden-Bell 1964; Larson 1975; Matteucci et al. 1998; Tantalo & 
Chiosi 2002), leads to the scenario of monolithic collapse, i.e. an earlier object formation 
and a following passive evolution. The clustering data can be used to discuss whether the 
merging processes were important at various redshifts or the galaxy number tends to be 
conserved during the evolution. These two opposite models predict a significantly different 
redshift evolution of the bias factor (Matarrese et al. 1997; Moscardini et al. 1998). 



1 Bascd on observations collected at the European Southern Observatory, Chile (ESO P66.A-0277 and ESO 
P67.A-0537), with the Arizona Steward Observatory and with National Telescope Galileo (TNG) during A03 
period. 
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The first attempt to measure the clustering of QSOs was made by Osmer (1981). Shaver 
(1984) was the first to detect QSO clustering on small scales using the Veron-Cetty & Veron 
(1984) catalog, a collection of inhomogeneous samples. A number of authors (Iovino & Shaver 
1988; Andreani & Cristiani 1992; Mo & Fang 1993; Shanks & Boyle 1994; Andreani et al. 
1994; Croom & Shanks 1996) have used complete and better defined QSO samples to measure 
their spatial distribution. At a mean redshift of z ~ 1.4, they generally detect a clustering 
signal at a typical significance level of ~ 3 — 4a, corresponding to a correlation length, r , 
similar to the value obtained for local galaxies: r ~ Qh^ 1 Mpc. However, there has been 
significant disagreement over the redshift evolution of QSO clustering, including claims for 
a decrease of r with redshift (Iovino & Shaver 1988), an increase of r with redshift (La 
Franca et al. 1998) and no change with redshift (Croom & Shanks 1996). Recently, Croom 
et al. (2001), using more than 10,000 objects taken from the preliminary data release catalog 
of 2dF QSO Redshift Survey (hereafter 2QZ), measured the evolution of QSO clustering as 
a function of redshift. Assuming an Einstein-de Sitter universe (Qm = 1-0 and Q\ = 0.0), 
they found no significant evolution for tq in comoving coordinates over the redshift range 
0.3 < z < 2.9, whereas for a model with Q M = 0.3 and Ha = 0.7 the clustering signal 
shows a marginal increase at high redshift. Here Qm and Q\ are the mass and cosmo logical 
constant density contribution, respectively, to the total density of the universe. 

The observed behavior of the QSO clustering can be explained within the linear theory 
and a typical bias model. The theoretical interpretation of the picture drawn by 2QZ is a 
result of the combination of many ingredients and their degeneracies: the bias factor, the 
ratio between the masses of black hole and dark matter halo, the life-time of QSOs, the 
efficiency and the mass accretion rate. 

To add new insights in the modeling and interpretation, one has to consider the con- 
straints from the luminosity function (LF) or/and to enlarge the redshift domain toward 
lower or higher redshifts. For these reasons, we have started a project, the AERQS, to find 
bright AGN in the local universe, removing present uncertainties about the properties of the 
local QSO population and setting the zero point for clustering evolution. For the general 
aims of the AERQS Survey and its detailed presentation [see Grazian et al. (2000, 2002)]. 

The goal of this paper is to analyze the clustering properties of a well defined large 
sample of bright QSOs at z < 0.3, and provide key information on the following issues: what 
is the typical mass of DMHs hosting AGN? What is the typical bias factor for AGN? What 
is the duty cycle for AGN activity? What is the typical efficiency of the central engine at the 
various redshifts? 

The plan of the paper is as follows. In § 2 we describe the data used in the statistical 
analysis. The various techniques used to investigate the clustering properties are presented 
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in § 3, while § 4 is devoted to a comparison with similar results obtained by previous surveys 
at low redshifts. To investigate the redshift evolution of the clustering, the spatial properties 
of QSOs in the local universe are compared in § 5 to the recent 2QZ results at intermediate 
redshifts for QSOs and to various estimates for normal and peculiar galaxies. In § 6 physically 
motivated models are used to link the galactic structures at high-z with the local AGN and 
galaxy population. In § 7 we discuss the clustering properties of QSOs in the light of these 
simple theoretical models. Finally § 8 gives some concluding remarks on the clustering of 
QSOs. 



2. The Data 

It is paradoxical that in the era of 2QZ and Sloan Digital Sky Survey (SDSS), with 
thousands of faint QSOs discovered up to the highest redshifts, there are still relatively few 
bright QSOs known at low redshift. One of the main reasons, as shown in previous papers 
(Grazian et al. 2000, 2002), is the rather low surface density of low- z and bright QSOs, of 
the order of few times 10~ 2 per deg 2 . This corresponds to a very small number of objects 
also in the case of the 750 deg 2 of the complete 2QZ and 1000 deg 2 of the SDSS (during 
commissioning phase). One more reason, not less important than the previous one, is that 
with the optical information only it is difficult to efficiently isolate bright QSOs from billions 
of stars in large areas. As a consequence, a survey based on different selection criteria is 
required. In Grazian et al. (2000, 2002) we have used the X-ray emission, a key feature of 
the AGN population. 

The AERQS is divided in three sub-samples, two in the northern hemisphere (USNO 
and GSC), described in Grazian et al. (2000) (hereafter Paper I), and the DSS sample in 
the southern hemisphere, described in Grazian et al. (2002) (Paper II). After a campaign of 
spectroscopic identifications at various telescopes, we have completed the sample, which is 
made up of 392 AGN with redshifts between 0.007 and 2.043. The redshift distributions, 
shown in Fig. 1, show a peak around z ~ 0.1 with an extended tail up to z = 0.4. Five AGN 
with 0.6 < z < 2.04 are possibly objects magnified by gravitational lensing effects. Table 1 
summarizes the basic properties of the three sub-samples. The area covered by the AERQS 
Survey consists of ~ 14,000 deg 2 at the high Galactic latitudes (|& ga i| > 30°). The mean 
values for the completeness and efficiency are 65.7% and 52.3%, respectively. 
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Fig. 1. — The redshift distributions for the three separated sub-samples (DSS, USNO and 
GSC) and for the total sample (AERQS, bottom-right panel). Five AGN with higher red- 
shifts (in the range 0.6 < z < 2.04) are not plotted here. 
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3. Measuring the Clustering in the AERQS 

The simplest way to analyze the clustering properties of a homogeneous and complete 
sample of QSOs is to compute the TPCF, £(r), in the redshift space. We choose to calculate 
£(r) for two representative cosmological models: (Q M ,Q\) = (1.0,0.0) and (0.3,0.7). We 
will call these cosmological models Einstein-de Sitter (hereafter EdS) and A, respectively. 

To compute £(r) we have used the minimum variance estimator suggested by Landy & 
Szalay (1993): 

_ QQ(r) - 2QR(r) + RR(r) 
Urj RR(r) ' 1 ' 

where QQ, QR and RR are the number of QSO-QSO, QSO-random and random-random 
pairs with a separation r± Ar, respectively. Here r is the comoving distance of two QSOs in 
the redshift space. We compute the TPCF in bins of Ar = 5/i -1 Mpc, where h is the Hubble 
constant, in units of 100 km s _1 Mpc -1 . The adopted values for the Hubble constant are 
h = 0.5 for the EdS model, and h = 0.65 for A. We generate 100 random samples and we 
use the mean values of QR(r) and RR(r) for the estimator. 

The correct generation of the random objects is in general the most critical aspect in the 
clustering analysis. This problem becomes fundamental in the case of a flux-limited sample, 
the AERQS. The area covered by our survey is not homogeneously distributed in the sky, 
due to the selection criteria adopted and the variable Galactic extinction. Consequently the 
"true" apparent magnitude limit of our survey is variable. Fig. 2 shows the effective area 
covered by the AERQS survey, limited by an exposure time t exp > 300 sec in the RASS-BSC 
(Voges et al. 1999) and at high Galactic latitudes (|6 ga i| > 30°). 

The results on clustering reported in this paper are derived by scrambling the redshifts, 
and the right ascension (RA) and declination (DEC) coordinates for the total (AERQS) 
sample. The random RA and DEC are derived from Fig. 2, while the redshifts are randomly 
extracted from the observed (Gaussian smoothed) redshift distributions (Fig. 1). 

In order to check the robustness of the results, we have carried out a more complex 
generation of random QSOs, which ensured the uniformity of the "synthetic" samples. The 
angular positions were chosen again randomly from the map in Fig. 2. Then, for each 
object we generated random values for redshift and absolute magnitude reproducing the LF 
estimated by La Franca & Cristiani (1997) and Grazian et al. (2000) in the redshift range 
0.04 < z < 2.2. In particular for $(Mb,z) we adopt a double power-law relation evolving 
accordingly to a Luminosity Dependent Luminosity Evolution (LDLE) model: 
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Fig. 2. — The black area shows the regions covered by the AERQS All Sky Survey after 
applying the selection criteria described in Paper I and Paper II. The projection is done here 
in RAcos(DEC) vs. DEC. 
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HM B ,z) 



10 0.4[M s -M*^)](a+l) + 10 0.4[M s -M*( 2 )](/3+l) 



(2) 



where 



M* B (z) = M* B (z = 2) - 2.5Hog[(l + z)/3] 



(3) 



and 



k = h + hi [M B - M B {z)]e~ z/0A if M B < M* B (z) 

k = h if M B > M* B {z) . 



The parameters a and (3 correspond to the faint-end and bright-end slopes of the optical 
LF, respectively, and M B (z = 2) is the magnitude of the break in the double power-law shape 
of the LF at z — 2. The actual values adopted in the LDLE parameterization, reported in 
Tab. 2, are derived by a fit to the observed LF. Extinction by Galactic dust is taken into 
account using the reddening E(B — V) as a function of position, calculated by Schlegel et 



This approach, though computationally expensive, avoids biases in the generation of 
random samples of QSOs. It reproduces the observed LF and the distribution of redshifts 
and apparent magnitudes. Fig. 3 shows the observed and random generated QSOs in the 
(z, Mb) space in the case of the EdS model. 

The results on the clustering of QSOs derived with this particular approach are con- 
sistent with the ones obtained with the scrambling of the redshifts. In the following all the 
computations will be carried out with the latter method. 

First, we calculate the TPCF integrated over a sphere, £(r), as a function of the sphere 
radius r, for the three sub-samples separately (DSS, USNO, GSC) and the total sample 
(AERQS). Fig. 4 reports the results for the EdS universe, while Fig. 5 refers to a A 
universe. The error bars in Fig. 4 and 5 represent the la interval for £(r) and are obtained 
by assuming a Poisson distribution (Gehrels 1986). 

In order to investigate the possible presence of a spurious clustering signal at large 
scales, we have computed the angular TPCF binned in intervals of 3 degrees (corresponding 
to ~ U.Qh- 1 Mpc comoving). Fig. 6 shows the absence of any significant bias on the large 
scales sampled by AERQS, up to 150 hr l Mpc. 

The signal shown in Fig. 4 and 5 for separations smaller than 15/i -1 Mpc is due to 25 and 
28 QSO pairs for the EdS and A models, respectively. For a completely random distribution, 
the expected number of pairs is 12 for EdS and 14 for the A model. Considering separations 



al. (1998). 
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Table 1. A summary of the AERQS Survey: the three sub-samples. 



Name DEC limit magnitude Area Nagn Redshift Completeness 



DSS -90<5<0 12.60 < B < 15.13 5660 111 0.012 < z < 0.680 0.63 
USNO < 5 < +90 13.50 < R < 15.40 8164 209 0.034 < z < 2.043 0.68 
GSC < 5 < +90 12.50 < V < 14.50 8164 72 0.007 < z < 0.573 0.63 



Note. — The reported area (in deg 2 ) is the fraction of the northern and southern hemispheres 
with |6 ga i| > 30° and exposure time of the ROSAT All Sky Survey t exp > 300 sec (as described 
in Paper I and Paper II). 



Table 2. The parameters used for the LF of QSOs. 



Model 


$* 


M* B (z = 2) 


a 


P 


h 


k 2 


EdS 
A 


9.8 
5.0 


-26.3 
-26.7 


-1.45 
-1.45 


-3.76 
-3.76 


3.33 
3.33 


0.37 
0.30 



Note. — $* is in units of 10 7 mag 1 Mpc 3 . 
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Fig. 3. — The redshift vs. magnitude Mb distribution for the QSOs of AERQS. The observed 
AGN sample (filled circles) is compared with the random generated sample (small dots) in 
the (z, Mb) space. The density of random points is 100 times larger than the observed one. 
Results are shown for the EdS model. 
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Fig. 4. — The integrated correlation function £(r), as a function of the sphere radius r, for 
an EdS model. Different panels refer to the integrated correlation function (and la error 
bars) for the DSS, GSC, USNO sub-samples, and for the AERQS total sample. £(r) is 
integrated over spheres of increasing radii. Consequently the error bars, which are shown 
only for reference, are not independent. At large scales (> 50/i _1 Mpc) the integrated TPCF 
is consistent with zero, showing the absence of large-scale gradients in the data. 
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Fig. 5. — As Fig. 4, but for the A model. 
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smaller than 20h 1 Mpc the observed pairs are 36 and 38, to be compared with 27 and 26 
random pairs expected. The clustering signal is therefore detected at a 3 — 4<r level. 

The differential TPCF for the complete AERQS sample is shown in Fig. 7, both for 
EdS (upper panel) and A models (lower panel). The results have been fitted by adopting a 
power-law relation 

« r) (4) 

The best-fit parameters can be obtained by using a maximum likelihood estimator 
(MLE) based on Poisson statistics and unbinned data (Croft et al. 1997). Unlike the usual 
X 2 -minimization, this method allows to avoid the uncertainties due to the bin size (see above), 
the position of the bin centers and the bin scale (linear or logarithmic). 

To build the estimator, it is necessary to estimate the predicted probability distribution 
of quasar pairs, given a choice for the correlation length r and the slope 7. The small number 
of pairs observed at small scales makes a reliable determination of the slope 7 particularly 
difficult. Therefore, we have used fixed values for the slope 7, adopting those obtained by 
Croom et al. (2001) for the 2dF catalog, namely 7 = 1.58 and 7 = 1.56 for EdS and A models, 
respectively. In this way, the comparison with the TPCF at higher redshifts obtained from 
the 2dF data is equivalent both in term of r and £. 

By using all the distances between the quasar-random pairs, we can compute the number 
of pairs g(r)dr in arbitrarily small bins dr and use it to predict the mean number of quasar- 
quasar pairs h(r)dr in that interval as 

N c — 1 

h(r)dr = ^ N i 1 + ^( r )]9(r)dr , (5) 

where the correlation function £ is modeled with a power-law as in Eq.(4) 2 . In this way, it is 
possible to use all the distances between the N p quasar-quasar pairs data to build a likelihood. 
In particular, the likelihood function C is defined as the product of the probabilities of having 
exactly one pair at each of the intervals dr occupied by the quasar-quasar pairs data and the 
probability of having no pairs in all other intervals. Assuming a Poisson distribution, one 
finds 

N p 

£ = Y\_ ex P[~ h(r)dr]h(r)dr J^[exp[— h(r)dr] , (6) 



2 Actually the previous equation holds only for the Davis & Peebles (1983) estimator (the original formu- 
lation for the TPCF, £(r) = — 1), but, since the results obtained using different estimators are similar, 
we can safely apply it here 
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where j runs over all the intervals dr where there are no pairs. It is convenient to define the 
usual quantity S = — 2 In C, which can be written, once we retain only the terms depending 
on the model parameter r , as 



The integral in the previous equation is computed over the range of scales where the 
fit is made. The minimum scale is set by the smallest scale at which we find QSO pairs 
(?"min = Mpc), while for the maximum scale we adopt r max = 30/i _1 Mpc. The latter 
choice is made to avoid possible biases from large angular scales, where the signal is weak. 

By minimizing S one can obtain the best-fitting parameter r . The confidence level is 
defined by computing the increase AS with respect to the minimum value of S. In particular, 
assuming that AS is distributed as a x 2 with one degree of freedom, AS* = 1 corresponds to 
68.3 per cent confidence level. It should be noted that by assuming a Poisson distribution 
the method considers all pairs as independent, neglecting their clustering. Consequently the 
resulting error bars can be underestimated [see the discussion by (Croft et al. 1997)]. 

In Fig. 7 the lines represent the la confidence region computed with the MLE method 
previously described, varying only the correlation length tq. We find r$ = 8.49^2 05^ _1 
Mpc for the EdS model (with 7 = 1.58) and r = 8.64^2Qg/i _1 Mpc for the A model (with 
7 = 1.56). The quoted errors on r are based on the assumption of a fixed slope. It is well 
known that the errors on r and 7 are correlated. Fixing the slopes to 7 = 1.58 and 1.56 
allows us to derive the confidence levels for the integrated TPCF £ which can be consistently 
compared with 2QZ results. 

It can be useful to present the previous results in a non-parametric form, specified by 
the clustering amplitude within a given comoving radius, rather than as a scale length which 
depends on a power-law fit to £(r). This is generally represented by the correlation function 
integrated over a sphere of a given radius in redshift-space r max , 



This is the same quantity we plotted in Figs. 4 and 5 for varying r max . Different authors have 
chosen a variety of values for r max , e.g. lOh^ 1 Mpc (Shanks & Boyle 1994; Croom & Shanks 
1996), 15/1- 1 Mpc (La Franca et al. 1998), or 20/i" 1 Mpc (Croom et al. 2001). In general, 
the larger the scale on which the clustering is measured, the easier the comparison with the 
linear theory of the structure evolution. Since in the following sections we will compare our 
results with those obtained for the 2QZ by Croom et al. (2001), we prefer to quote clustering 






max 



(8) 



- 16 - 



1 1 1 1 1 1 1 1 1 

- 0.02 < z < 0.22 

n„=i o n =o o 

- M A 

- DSS 


i i | i i i | i i i 


i i i i i i i i i i i i i i i i i i i 

- 0.04 < z < 0.22 

n =i o n =o o 

- M A 

1 USNO 1 






\ ^\\\^\\\\\\\\\# _ 


_l 1 1 1 1 1 1 

- 0.02<z<0.22 

n =1.0 n =o.o 

M A 

Lgsc 


i — i — | — i — i — i — | — i — \—v 


_l 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

- 0.02<z<0.22 

n u =i.o n =o.o 

1 AERQS J 


r i i 1 i i i 1 i 


i i i i i i i i i i 


, , , i , , , i , , , i , , , i , , , 



20 40 60 80 20 40 60 80 

©(degree) ©(degree) 



Fig. 6. — The differential angular TPCF binned in intervals of 3 degrees (corresponding to 
~ 12.6/2T 1 Mpc). Different panels refer to the correlation function (and la error bars) for 
the DSS, GSC, USNO sub-samples, and for the AERQS total sample. These results do not 
depend on the adopted cosmological model. The angular TPCF at small scales is consistent 
with zero because it is diluted over 10-20 hr x Mpc. 
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Fig. 7. — The (differential) two-point correlation function (with la error bars) for the AERQS 
sample in EdS (top panel) and A models (bottom panel). The solid lines show the la 
confidence region obtained by fitting the data with a power-law relation (with fixed slope) 
using the MLE approach described in the text. Only points at r < 30/i _1 Mpc are used to fit 
the TPCF, and they are indicated by the solid line. The dashed line indicates the extension 
of the TPCF relation to data that are not used for the fit. At small scales (log(r) < 1.8) a bin 
size of 5/i -1 Mpc has been adopted, while at larger scales (log(r) > 1.8) a bin sizes of lOh -1 
Mpc has been used. The bin size was chosen too avoid too large bin to bin fluctuations. 
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amplitudes within 20h~ 1 Mpc, a scale for which linearity is expected to better than a few 
per cent. Choosing a large radius also reduces the effects of small scale peculiar velocities 
and redshift measurement errors, which may well be a function of redshift. 

Table 3 summarizes the values of r , 7 and £(20) for the total sample, both for the 
EdS and A models. In the same table, we list the mean redshift (z) of the observed QSO 
sample. We also report the median value of the redshift of the QSO pairs computed within 
a sphere of 20/i^ 1 Mpc, z%. We find that it is systematically lower than the mean redshift of 
the sample. 

In our analysis, we do not take into account the velocity field of QSOs, the cone edge 
effect and the effect of statistical errors on QSO redshifts. Recent papers [see e.g. Croom 
et al. (2001)] suggest that the Poisson errors, due to the limited size of a sample, are more 
important than these effects. 

4. Comparison with other surveys 

It is instructive to compare the present results on the clustering of low- z AGN with 
that of other surveys, both at low- and high- redshift, in order to get information about the 
connection between various galactic structures and their evolution. To avoid problems with 
different assumptions on the values of the slope 7, we decided to compare the values of the 
integrated TPCF at 20 /i -1 Mpc, £(20). When not directly available in the original paper, 
£(20) has been computed by integrating the TPCF with the best fitting values of r and 7. 

4.1. Comparison with other local AGN surveys 

Using a low-redshift (z < 0.2) sample, Boyle & Mo (1993) measured the clustering 
properties of 183 AGN in the EMSS. They found evidence for a small value of the integrated 
TPCF, £ = 0.7 ± 0.6 (computed at lOh" 1 Mpc), corresponding to a correlation length of 
r = b.Otlih- 1 Mpc. The assumed slope for the TPCF is 7 = 1.8 and the resulting £(20) 
is 0.20 ± 0.17. Considering the uncertainties, this result is slightly lower than or consistent 
with our results. Moreover, since the Boyle & Mo sample is obtained by identifications of 
X-ray sources, it contains fainter 3 AGN than AERQS. As a consequence, a slightly smaller 
value of ro is expected for their sample, because the clustering strength is found to depend, 
weakly, on the absolute magnitude M B , as shown in Croom et al. (2002) and in Norberg et 



3 AGN in the EMSS are typically 5 times fainter than L* at z ~ 0.2, or 1.75 magnitude fainter than Mg. 
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al. (2002). 

Georgantopoulos & Shanks (1994) investigated the clustering properties of 192 Seyfert 
galaxies from the IRAS All Sky Survey. They claimed a 2 — 3a detection at 10 — 20h~ 1 
Mpc, corresponding to £(20) = 0.14 ± 0.15 at z — 0.05, similar to local late-type galaxies. 
This result is consistent with a model in which local QSOs randomly sample the galaxy 
distribution. 

Carrera et al. (1998) analyzed the clustering of 235 X-ray selected AGN with < z < 3, 
obtaining an integrated TPCF of 0.02 < £(20) < 0.25. The redshift range of this survey 
is particularly extended and the density of sources correspondingly low. The clustering 
detection is marginal, at 2a level only. Moreover there are only 33 AGN with z < 0.2 in this 
sample. 

Akylas et al. (2000) investigated the angular correlation function of 2096 sources selected 
from the RASS-BSC. They rejected known stars and other contaminants: a cross-correlation 
analysis with spectroscopic samples indicated that the majority of their sources are indeed 
AGN. They obtained a ~ 4<r detection of clustering. Using the Limber equation and assum- 
ing a source redshift distribution (not shown in their paper) with an estimated mean value of 
0.1, they derived £(20) = 0.35 ± 0.09. Stars, galaxy clusters or other spurious contaminants 
could affect their results. 

Mullis et al. (2001) derived the clustering properties of 217 AGN found in the North 
Ecliptic Pole (NEP) survey, a connected area of ~ 81 deg 2 covered by ROSAT observations. 
The sample spans the redshift interval < z < 3.889, with z = 0.408. A 3.8<x clustering 
detection was obtained, corresponding to an integrated TPCF of £(20) = 0.36 ± 0.15. This 
result confirms that X-ray selected AGN are spatially clustered in a manner similar to that 
of optically/UV selected AGN. 

Notice that Boyle & Mo (1993), Georgantopoulos & Shanks (1994), Carrera et al. (1998) 
and Akylas et al. (2000) used an EdS cosmology to compute the clustering properties of their 
samples, while Mullis et al. (2001) adopted a A model. 

Finally, it is interesting to compare the clustering properties of AGN and normal galaxies 
at low-2, using our results and recent results by Norberg et al. (2002). Our value for the 
AGN correlation strength (£(20) ~ 0.461 ± 0.237) appears slightly lower than the typical 
value for the brighter early-type galaxies (£(20) = 0.70^0^), or at most consistent, indicating 
that these two classes have not experienced a completely different evolutionary history, but 
could represent two distinct phases during the processes of formation and evolution of the 
same objects. This gives an additional support for models dealing with the joint evolution 
of QSOs and normal galaxies [e.g. see Haehnelt & Kauffmann (2000); Granato et al. (2001); 



-20- 



Franceschini et al. (2002) and references therein]. In particular, the fact that the correlation 
length of AGN at z ~ is consistent with that of ellipticals or SO in the local universe, 
reinforces the hypothesis that the QSO host galaxy should be old. 

4.2. Comparison with QSO clustering at high-redshift 

La Franca et al. (1998) investigated the evolution of QSO clustering using a sample 
of 388 QSOs with 0.3 < z < 2.2 over a connected area of 25 deg 2 down to B < 20.5 
magnitude. Evidence was found for an increase of the clustering with increasing redshift 
(£(20) = 0.22+°;^ at 0.3 < z < 1.4 and £(20) = 0.87^ at 1.4 < z < 2.2). This result 
does not support the idea of a single population model for QSOs. The general properties of 
the QSO population studied by La Franca et al. (1998) would arise naturally if QSOs are 
short-lived events (r ~ 10 6 — 10 7 yr) related to a characteristic halo mass of ~ 5 ■ 10 12 M Q . 

Croom et al. (2001) have used more than 10,000 QSOs taken from the preliminary data 
release catalog of 2QZ to measure the QSO clustering as a function of redshift. Their sample 
spans two connected areas for a total of 750 deg 2 at a limiting magnitude of 20.85 in the Bj 
band. The completely identified sample (not yet released) consists of nearly 22,500 QSOs 
in the redshift range 0.3 < z < 2.2. The results from the preliminary data release (to 
be considered with some caution), expressed in terms of the correlation function integrated 
inside spheres of 20/i -1 Mpc, £(20), are shown in Fig. 8, together with the estimates obtained 
at z ~ 0.1 for the AERQS sample. The discussion of the theoretical models shown by the 
different lines, will be given in the following Section. 

For an EdS universe (left panels), Croom et al. (2001) find that there is no significant 
evolution of the QSO clustering in comoving coordinates over the whole redshift range con- 
sidered. Assuming a A model (right panels), the clustering shows a marginal increase at 
high redshifts, with a minimum of £(r) near z ~ 0.5 — 1.0. Our data show a tendency to an 
increase of £(20) at low-z, both for the EdS and A models. This result supports the predic- 
tions based on simple theoretical models (see next section) and on numerical simulations by 
Bagla (1998). Notice that very recently this general trend for the clustering evolution has 
been also confirmed by the power spectrum analysis made by Outram et al. (2003) using the 
final version of the 2QZ catalog, containing 22,652 QSOs. 

The AERQS AGN catalog samples a part of the QSO luminosity function which is 
fainter than that sampled by the 2QZ. The mean absolute magnitudes of the total AERQS 
QSOs are M B = -23.49 and -22.99 for the EdS and ACDM, respectively. The 2QZ QSOs 
at 0.3 < z < 2.2 are brighter than local AGN, with Mj = -24.43 and -25.11 for EdS and A, 
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respectively. In our comparison, we do not take into account the dependence of the TPCF on 
the absolute magnitude of the sample, since the QSO population exhibits a strong luminosity 
evolution with redshift and Croom et al. (2002) have demonstrated that the dependence of 
the clustering on Mb is very weak. A correction of the luminosity dependence of the TPCF 
would increase the AERQS value, enhancing the redshift evolution of the clustering. 



5. Modeling the redshift evolution of QSO clustering 

By adding the AERQS value of local QSO correlation length to the 2QZ estimates at 
higher redshifts, we have now the complete picture of the QSO clustering properties up to 
z ~ 2.5, as summarized in Fig. 8. In the following, we introduce a model which can be used 
to interpret the observed evolution. 

In general, the theoretical understanding of how matter clustering grows via gravita- 
tional instability in an expanding universe is presently quite well developed, even if the 
number of ingredients required in the models is large. As a consequence, it is relatively 
straightforward to compute the correlation function of matter fluctuations, £ m , as a function 
of redshift, given a cosmological scenario [see e.g. Peacock & Dodds (1996); Smith et al. 
(2003)]. However, this does not lead directly to a prediction of QSO correlation properties 
because the details of the link between the distribution of active nuclei and the distribution 
of the mass are not fully understood. In principle, this relationship could be highly complex, 
non- linear and environment-dependent, making very difficult to obtain useful informations 
on the evolution of matter fluctuations from the AGN clustering. In this spirit, a relatively 
simple form of the local bias b is generally assumed. 

Matarrese et al. (1997; see also Moscardini et al. 1998; Hamana et al. 2001) developed 
an algorithm for describing the clustering on our past light-cone taking into account both 
the non-linear dynamics of the dark matter distribution and the redshift evolution of the 
bias factor. The final expression for the observed spatial correlation function £ b s in a given 
redshift interval Z is 

f z dzidzJJ(zi)Jf(z2) b cS (z 1 )b cS (z 2 )U(r, z) 
tobs(r) — , (9 J 

IJz dziNiz!)} 

where JJ"(z) = J\f(z)/r(z), M{z) is the actual redshift distribution of the catalog and r{z) 
describes the relation between comoving radial coordinate and redshift. Here z is a suit- 
ably defined intermediate redshift. The method has been extended to include the effects of 
redshift-space distortions using linear theory and the distant-observer approximation (Kaiser 
1987). 
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Fig. 8. — The evolution of the integrated TPCF, £(20). This plot compares the observational 
data, shown by points with la error bars, to the predicted redshift evolution of the clustering 
for QSO-conserving (top panels) and merging (bottom panels). The point at z = 0.1 is our 
AERQS result, while points at higher redshifts come from 2QZ analysis. Theoretical results 
for SCDM and ACDM are shown in the left and right panels, respectively. The dashed lines 
represent the clustering evolution for models with a given value of the correlation length ro 
at z — 0, as indicated in the plot. The solid line in the bottom- right panel refers to the 
predictions of a combined model, which assumes a merging phase (with logM min = 12.5) at 
high-2 and a following conserving phase at low-z (see the text for more details). 
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A fundamental role in the previous equation is played by the effective bias 6 e fr- In fact 
the final aim of models dealing with clustering is to determine the behavior of the bias factor, 
once a given theoretical picture is assumed. In practice the effective bias can be expressed 
as a weighted average of the 'monochromatic' bias factor b(M, z) of objects with some given 
intrinsic property M (like mass, luminosity, etc): 



J M 

where Af(z, M) is the number of objects actually present in the catalog with redshift within 
dz of z and property within rflnM of InM, whose integral over InM is Af(z). 

In most fashionable models of structure formation, the growth of large-scale features 
happens because of the hierarchical merging of sub-units. Since the development of the 
clustering hierarchy is driven by gravity, the most important aspects to be understood are 
the properties of dark halos rather than the QSOs residing in them. Following Mo & White 
(1996), it is possible to calculate the bias parameter b(M,z) for halos of mass M and 'for- 
mation redshift' Zf observed at redshift z < Zf in a given cosmo logical model as 



where a 2 M is the linear variance averaged over the scale corresponding to the mass M, extrap- 
olated to the present time [z — 0); 8 C is the critical linear over-density for spherical collapse; 
D + is the growing factor, depending on the cosmo logical parameters Qm and The dis- 
tribution in redshift and mass n(z, M) for the dark halos can be estimated using the Press 
& Schechter (1974) formalism; in particular in the following analysis we adopt the relation 
found by Sheth & Tormen (1999). In the standard treatment of hierarchical clustering, all 
the halos that exist at a given stage merge immediately to form higher mass halos, so that 
in practice at each time the only existing halos at all are those which just formed at that 
time (i.e. Zf — z). If one identifies quasars with their hosting halos, then the merging rate 
is automatically assumed to be much faster than the cosmological expansion rate. This is at 
the basis of what Matarrese et al. (1997) and Moscardini et al. (1998) called merging model. 
Of course this instantaneous-merging assumption is physically unrealistic and is related to 
the fact that we use a continuous mass variable, while the aggregates of matter that form are 
discrete. Assuming a monotone relation between the mass and the observational quantity 
defining the limits of a given survey, the effective bias can be estimated by considering that 
the observed objects represent all halos exceeding a certain cutoff mass M min at any partic- 
ular redshift. In this way, by modeling the linear bias at redshift z for halos of mass M as 
in equation (11) and by weighting it with the theoretical mass-function n(z, M), which can 
be self-consistently calculated using the Sheth & Tormen (1999) relation, the behavior of 




(10) 




(11) 
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b e g(z) is obtained. The parameter M min can be regarded as a free parameter or alternatively 
fixed in order to obtain given values of the correlation length r at z = (see later). 

An alternative picture of biasing can be built by imagining that quasar formation occurs 
at a relatively well-defined redshift Zf. Actually there are no changes if one assume that there 
is some spread in the distribution of Zf. If this is the case, one can further imagine that 
quasars, which are born at a given epoch Zf, might well be imprinted with a particular value 
of b(M, Zf) as long as the formation event is relatively local. If quasars are biased by birth in 
this way, then they will not continue with the same biasing factor for all time, but will tend 
to be dragged around by the surrounding density fluctuations, which are perhaps populated 
by objects with a different bias parameter. In this case, the evolution of the bias factor can 
be obtained from (Fry 1996): 

b(z) = l + (b f -l)2±&, z<z f , (12) 

where bf is the bias at the formation redshift Zf. Notice that b(z) approaches unity with 
time, provided that the universe does not become dominated by curvature or vacuum in the 
meantime (Catelan et al. 1998). This model is called conserving model [see Matarrese et 
al. (1997) and Moscardini et al. (1998)] or, alternatively, test particle model. Again, it is 
difficult to motivate this model in detail because it is hard to believe that all galaxies survive 
intact from their birth to the present epoch, but at least it gives a plausible indication of the 
direction in which one expects b to evolve if the timescale for quasar formation is relatively 
short and the timescale under which merging or disruption occurs is relatively long. 

Notice that the merging model (rapid merging) and conserving model (no merging) can 
be regarded as two extreme pictures of how structure formation might proceed. In between 
these two extremes, one can imagine more general scenarios in which quasars neither survive 
forever nor merge instantaneously. The price for this greater generality is that one would 
require additional parameters to be introduced in the models (see the discussion at the end 
of the next section). 



5.1. Results 

In the following analysis we will present the results for two different cosmological models. 
Both models assume a CDM power spectrum (Bardeen et al. 1986), with spectral index n = 1 
and shape parameter T = 0.2. The power spectrum normalization (expressed in terms of 
as, i.e. the r.m.s. fluctuation amplitude in a sphere of 8k" 1 Mpc) is chosen to be consistent 
with very recent estimates obtained from the cluster abundance analysis [e.g. Reiprich & 
Bohringer (2002); Viana et al. (2002); Seljak (2002)]. The two considered models are: 
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• a "standard" CDM Einstein-de Sitter model with a 8 = 0.5 (hereafter SCDM); 

• a flat CDM universe with (n M ,Sl A ) = (0.3,0.7) with a 8 = 0.8 (hereafter AC DM). 

In Fig. 9 we show the redshift evolution of the bias for different values (indicated in 
the plot in units of /i~ 1 M ) of the minimum mass M min of the halos hosting QSOs. The 
theoretical predictions are compared to the observational results, shown by the points with 
la error bars. The values at z — 0.1 represent the bias parameter derived by our analysis 
of the AERQS. They are obtained by dividing the measured integrated TPCF for QSOs, 
£(20) by the theoretically predicted autocorrelation function of the underlying matter £ m (20): 
b 2 = £(20)/£ m (20). We obtain b = 1.75 ± 0.51 and b = 1.37 ± 0.35 for SCDM and ACDM, 
respectively. The data for higher redshifts come from the analysis of the 2QZ survey (Croom 
et al. 2001). 

A first comparison shows that the values for AERQS are consistent with the values at 
z ~ 0.7 for the 2QZ, implying the absence of a significant evolution of bias at low redshifts. 
As already noticed by Croom et al. (2001), the trend at higher redshifts for bias appears in 
general to depend on the cosmological models: for ACDM model, the observed b is always 
an increasing function of redshift, while in the SCDM case the value of b is almost constant 
for z > 1.5. 

More interesting is the comparison of the observed b with the theoretical predictions 
obtained assuming different M min . For the ACDM model the AERQS value corresponds 
to logM min = 12. 7+^7 (la error bars), and the observed trend is consistent with the 
bias evolution expected for dark halos with a minimum mass almost constant in redshift 
(logM min ~ 12 — 12.5). On the contrary, for the SCDM model it is impossible to repro- 
duce the bias factor using a constant minimum mass: while the value for AERQS suggests 
logM min = 12.5 ± 0.7 (always la error bars), the bias factor corresponds to halos with 
logM min ~ 11.5 at 0.5 < z < 1.5 and logM min < 11.5 at z > 2. 

In Fig. 8 we show the predictions for the redshift evolution of the TPCF integrated over 
20h~ 1 Mpc computed adopting the Q SO- conserving (upper panels) and merging (bottom 
panels) models, described above. The points (with la error bars) refer to the observational 
estimates, again from AERQS at z — 0.1 and from 2QZ at higher redshifts. The dashed lines 
represent the results obtained for models built to have given values of the QSO correlation 
length r at z = 0. In particular in the case of the QSO-conserving model, we show results 
for r (z = 0) = 6, 8, 10 h" 1 Mpc for both models. In the case of the merging model, we show 
the results for r (z = 0) = 4, 6, 8/i _1 Mpc, corresponding to a minimum dark matter halo 
mass of 3.4 x 10 n ,4.6 x 10 12 , 1.7 x lO 13 /!- 1 ^/©, for SCDM, and for r (z = 0) = 6,8, lOh' 1 
Mpc, corresponding to M min = 2.5 x 10 12 , 1.2 x 10 13 /i~\2.9 x 10 13 h^M & for ACDM. 
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Fig. 9. — The redshift evolution of the QSO bias factor. The points with la error bars 
represent the observational estimates: the point at z = 0.1 comes from this analysis of the 
AERQS catalog, while the remaining data are from 2QZ. The lines show the evolution of 
the bias obtained assuming different values of the mass M min (in units of /i _1 M ). 
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From the figure, it is evident that in the case of SCDM the QSO-conserving model is 
more or less able to reproduce the clustering evolution over the whole redshift interval, once 
a local value of r ~ Ih^ 1 Mpc is used. The situation is quite different for ACDM, for which 
the high clustering observed at z > 2 is not compatible with any trend predicted by the 
QSO-conserving model: only for z < 1 the decrease of £(20) follows the model expectations 
corresponding to a local value of tq{z = 0) ~ 7/i -1 Mpc. 

The comparison of the QSO clustering with the predictions of the merging model shows 
only marginal agreement on the whole interval 0.0 < z < 2.5. In particular, for SCDM the 
observational data are close to the predictions of the model corresponding to a low value of 
the local clustering, ro(z = 0) ~ 4/i _1 Mpc, while for ACDM the better agreement is for 
models corresponding to r (z = 0) ~ 7/i~ 1 Mpc, but with large deviations. 

As already said, these simple schemes do not exhaust all the possible scenarios through 
which QSOs might have formed and evolved. For example, it is quite possible that merging 
could play a different role at different redshifts. Present-day AGN, for example, have clearly 
not just formed at the present epoch since their observational properties suggest a lack of 
mergers in the recent past. On the other hand, it is plausible that QSOs at much higher 
redshifts, say z > 2, are undergoing merging on the same timescale as the parent halos. This 
suggests the possible applicability of a model where rapid merging works at high redshifts, 
but it ceases to dominate at lower redshifts and the bias then evolves by equation (12) 
until now. In this context it is interesting to note that, while bf is a free parameter in 
equation (12), it is actually predicted, once the appropriate minimum mass is specified. In 
fact it corresponds to the bias at the redshift Zf when objects stop merging. This model has 
been introduced by Moscardini et al. (1998), where the resulting relations for the redshift 
evolution of the bias factor are given. The complete application of this combined model 
to the present data on quasar clustering is quite difficult because of the size of error bars. 
Only as an example, for ACDM we compute the predicted £(20) by assuming a merging 
model with logM min = 12.5, followed by a conserving phase. The result, shown as solid 
line in the bottom- right panel of Fig. 8, is in rough agreement with the observational data, 
indicating that the redshift Zf, where the transition between the two different regimes occurs, 
is located at approximately Zf = 0.8. Of course, a validation of this model requires more 
robust estimates of the QSO clustering properties. 
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6. Discussion 



6.1. Clustering in the local universe 



As discussed in the previous paragraphs, there are empirical and theoretical evidences 
that QSOs are biased with respect to the matter distribution. The results of Boyle & Mo 
(1993), Georgantopoulos & Shanks (1994), Carrera et al. (1998), Akylas et al. (2000) and 
Mullis et al. (2001) are consistent, within the uncertainties, with the present results. In 
general we find that, at low-z, AGN have a correlation length of ~ 8/i" 1 Mpc. This value is 
quite similar to the correlation length of ellipticals, EROs or RGs at z ~ and higher than 
that of spirals or late-type galaxies. Assuming the hierarchical clustering paradigm and the 
SCDM model, a local correlation length of ~ 8.5±2/i _1 Mpc corresponds to a population of 
DMHs with mass larger than ~ 10 12 - 5 /i~ 1 M o (10 118 - lO^/i^M© at a la confidence level) 
and a bias parameter of b ~ 1.7. The space density of DMHs more massive than this limit is 
Pdmh ~ 3.7M0~ 3 /i 3 Mpc" 3 (0.63-19.64- 10~ 3 /i 3 Mpc" 3 at kr), as obtained by applying the 
Press-Schechter formalism (e.g. Sheth & Tormen 1999). The space density of bright AGN 
in the local universe is Pagn ~ 5.7- 10~ 7 /i 3 Mpc" 3 , as inferred by Grazian et al. (2000) using 
a sub-sample of the AERQS with limiting magnitude Mb = —22.5. We can thus obtain a 
rough estimate of the duty-cycle of local AGN, tagn using the simple relation 



where th is the Hubble time 4 . The duty-cycle of AGN at z ~ 0.1 turns out to be tagn ~ 
1.7-10 6 yr (the la confidence region is 3.3-10 5 — 1.0-10 7 yr). This result is in good agreement 
with the one Kauffmann & Haehnelt (2002) obtained by comparing their model to the 2QZ 
data, and only marginally consistent with the value of tqso ~ 10 7 yr obtained by Steidel et 
al. (2002) for QSOs at z ~ 3. It is worth noting that both these papers adopt a A model 
and h = 0.7. 

If we assume the relation found by Ferrarese (2002), namely 



it is possible to infer the mass of active BHs M B h at z ~ 0.1. With the values of the AGN 
DMH mass estimated in our sample, log/iM min /M Q = 12.5 ± —0.7, we can obtain a rough 
estimate for M BH : 6.7- lO 7 /^ 1 M (4.7- 10 6 - 9.5 • lO 8 /i _1 M at la). Assuming an absolute 



4 th is computed at z = 0.1 and corresponds to 11.3 -10 9 yr and 13.1 -10 9 yr for the EdS and A cosmological 
models, respectively 
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magnitude of M B = —22.5 that corresponds 5 to a bolometric luminosity L Bo i = 1O 12 L , one 
can infer a typical ratio L/M for a local AGN of 1.5 • 1O 4 L /M (1.1 • 10 3 - 2.1 • 1O 5 L /M ). 
For comparison, the Eddington value is L/M = 3.5 • 1O 4 L /M . In this case an efficiency of 
7] = L/L Edd ~ 0.4 (0.03 - 6.1) is derived. 

For the ACDM model, the correlation length is r = 8.6±2.0/i _1 Mpc, which correspond 
to a DMH mass of the order logM min = 12.7t° 8 r Following the same approach carried out 
for the SCDM model, with an AGN density of p A GN ~ 4.9 • 10 _7 /i 3 Mpc~ 3 and p DMU ~ 
8.12 • 10- 4 h 3 Mpc' 3 (1.28 • 10~ 4 - 4.31 • 10~ 3 /i 3 Mpc~ 3 at la) we derive a slightly longer 
duty cycle of t AG n ~ 7.9 • 10 6 yr (1.5 • 10 6 - 5.0 • 10 7 yr) and an efficiency of 4.8 • 1O 3 L /M 
(3.3 • 10 2 - 1.0 • 1O 5 L /M ), corresponding to rj ~ 0.14 (0.01 - 2.8, at la). Due to the 
large error bars on the TPCF, the constraints on the efficiency 77 are not stringent, but give 
nevertheless an indication of the mean value of the Eddington ratio for local AGN. 

AGN in the AERQS sample seem thus to accrete in a sub-Eddington regime, lower 
than the nearly- or super-Eddington accretion generally assumed for QSOs at high redshifts. 
High-2 QSOs are thought to have relatively small masses, so their extreme luminosities 
point to a high Eddington ratio (1-10 of the standard value). The direct determination 
of the Eddington ratio for QSOs at high-z is complicated by a number of difficulties, as 
discussed by Woo & Urry (2002), who suggest that the true value at z > 1 is uncertain and 
dominated by selection effects. At z ~ 0.1 they obtain a value LBoi/^Edd ~ 0.1, consistent 
with our result. Bechtold et al. (2003), using Chandra observations of high-redshift QSOs, 
estimated the BH mass and the Eddington ratio at 3.7 < z < 6.28 and compared it with 
the value for local AGN. At high-z QSOs possess masses of the order of 10 10 /i _1 M and are 
growing at a mass accretion rate of 0.1 mEdd- At low- z their results are comparable with our 
values, with M BH between 10 8 and lO 9 ^" 1 M and an Eddington ratio 77 between 10~ 2 and 
10 _1 . From the point of view of the theoretical modeling, Ciotti et al. (2003) were able to 
reproduce the QSO LF the mass function of local BHs with an Eddington ratio constant and 
equal to 10" 1 in the redshift range < z < 4. Haiman & Menou (2000) used an Eddington 
ratio decreasing from z = 4 to z = with a typical value of 10~ 2 — 10~ 3 at z ~ 0. 



6.2. Interpreting QSO and galaxy clustering 

In a sense, the general picture emerging from the observational data is that the clustering 
evolution for galaxies is similar to the QSO one: at low-z, the correlation length is decreasing 
from the local value reaching a minimum at z ~ 1, then it increases till z ~ 3—4. As discussed 



5 assuming a ratio Lboi/^b = 10. 
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by Arnouts et al. (1999, 2002), the term "evolution" has not to be considered literally. Given 
a survey defined by its characteristic limiting magnitude and surface brightness, the galaxies 
observed at high-z typically have higher luminosities. Therefore, the intrinsic differences of 
the galaxy properties at different z can mimic an evolution, i.e. the evolution measured in a 
flux-limited survey is not only due to the evolution of a unique population but can be due to 
a change of the observed population. The picture emerging from QSO clustering points in 
the same direction: QSOs are not part of a unique population, due to their short duty-cycle, 
and are intrinsically related with galaxy evolution. 

At this stage, it is important to discuss how the clustering depends on absolute magni- 
tude. Moreover, we should consider how the bias factor changes when the catalog selection 
effects are considered, i.e. when the theoretical quantities, as the mass M, are substituted 
by the observational ones, such as the luminosity L. Chosen one of the previously described 
models, one will end up with the quantity b(M, z) to be understood as 'the bias that ob- 
jects of mass M have at redshift z\ The effective bias at that redshift can be written more 
precisely as (see also Martini & Weinberg 2001) 



where N(z) = f dlnL $ obs (L) and <3> b s (£) is the observed luminosity function of the catalog, 
i.e. the intrinsic luminosity function multiplied by the catalog selection function, which will 
typically involve a cut in apparent magnitude, whatever wave-band is being used. 

Croom et al. (2002) observed a weak trend of the clustering strength with the magnitude, 
brighter objects being more clustered. Such a behavior can be understood by focusing on the 
results obtained in the previous paragraph. The efficiency and the accretion rate for local 
AGN determine the relation between mass and luminosity, and how they evolve with redshift. 
As a simple consequence, for QSOs at different epochs, luminosity does not necessarily trace 
mass. The dependence of clustering strength on the luminosity is thus weaker than the one 
expected in theoretical models assuming a fixed M/L ratio. 

EROs and RGs at z ~ 1 have higher correlation amplitude than QSOs at the same 
redshift; indeed they are following slow or passive evolution. Probably at z ~ they will 
plausibly become the brightest and most massive galaxies inside the clusters. QSOs, instead, 
show a different evolution for the TPCF: their behavior is consistent with that of a typical 
merging model at high-z and a passive evolution or object-conserving at \ow-z. 

Kauffmann & Haehnelt (2002) explored theoretically the possibility of using the cross- 
correlation between QSOs and galaxies, ^QGah to obtain new information on the masses of 
DMHs hosting QSOs. They used a semi-analytical model in which super-massive BHs are 




(15) 
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formed and fueled during major mergers. The resulting DMH masses can be in principle 
used to estimate the typical QSO life-time. In current redshift surveys, like the 2dFGRS 
or SDSS, these measurements will constrain the life-times of \ow-z QSOs more accurately 
than QSO auto-correlation function, because galaxies have much higher a space density than 
QSOs. As a result, ^,QGai can yield information about the processes responsible for fueling 
SMBHs. 



7. Conclusions 

The Asiago-ESO/RASS QSO survey (AERQS), an all-sky complete sample of 392 spec- 
troscopically identified objects (B < 15) at z < 0.3, has been used to carry out an extended 
statistical analysis of the clustering properties of local QSOs. 

The AERQS makes it possible to remove present uncertainties about the properties of 
the local QSO population and fix an important zero point for the clustering evolution and 
its theoretical modeling. With such a data-set, the evolutionary pattern of QSOs between 
the present epoch and the highest redshifts is tied down. 

On the basis of the (integrated and differential) two-point correlation functions, we have 
detected a 3— 4<r clustering signal, corresponding to a correlation length r = 8.6±2.0/i _1 Mpc 
and a bias factor b = 1.37±0.35 in a ACDM model. A similar value of ro, but corresponding 
to b — 1.75 ± 0.51, is obtained for an Einstein-de Sitter model, confirming previous analysis 
(Boyle & Mo 1993; Georgantopoulos & Shanks 1994; Carrera et al. 1998; Akylas et al. 2000; 
Mullis et al. 2001). These results shows that low-redshift QSOs are clustered in a similar way 
to radio galaxies, EROs and early-type galaxies, while the comparison with recent results 
from the 2QZ at higher redshifts shows that the correlation function of QSOs is constant in 
redshift or marginally increasing toward low redshifts. 

This behavior can be interpreted with physically motivated models, taking into account 
the non-linear dynamics of the dark matter distribution, the redshift evolution of the bias fac- 
tor, the past light-cone and redshift-space distortion effects. The application of these models 
allows us to derive constraints on the typical mass of the dark matter halos hosting QSOs: 
we have found logMoMH — 12.71q7 (the- mass is units of /i -1 M ), almost independently of 
the cosmological model. Using the abundance of dark matter halos with this minimum mass 
and assuming the relation found by Ferrarese (2002) between the masses of dark matter 
halos and active black holes, from the clustering data we can directly infer an estimate for 
the mass of the central active black holes and for their life-time, M B h ~ 2.1 • 1O 8 /i _1 M 
(1.0 • 10 7 - 2.9 • lOPh^Mo) and r AGN ~ 7.9 • 10 6 yr (1.5 • 10 6 - 5.0 • 10 7 yr), respectively. This 
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means that local AGN seem to accrete in a sub-Eddington regime. All these values have 
been obtained for a ACDM model; slightly shorter duty cycles are derived for an Einstein-de 
Sitter model. The time-life of z ~ 3 QSOs is ~ 10 7 yr, measured by Steidel et al. (2002). 
This could be a first indication that QSOs at all epochs have a similar life-time, which does 
not depend strongly on the Hubble time. 

Observational data suggest that most nearby galaxies contain central super-massive 
black holes, supporting the idea that most galaxies pass through a QSO/AGN phase. How- 
ever, the different clustering properties, together with the the short lifetimes derived for local 
AGN, suggest that this phase picks out a particular time in the evolution of galaxies, e.g. 
epochs of major star formation, interactions or merging. In this way the study of the QSO 
clustering evolution using extended catalogs helps us to distinguish between a number of 
possible QSO formation mechanisms. 
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Table 3. A summary of the clustering properties of the AERQS Sample. 
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(1.0,0.0) 


8.49 


6.44- 


-10.46 


1.58 


0.089 0.063 


0.368 


0.151- 


-0.585 


1.75 ±0.51 


(0.3,0.7) 


8.64 


6.56- 


-10.64 


1.56 


0.088 0.062 


0.461 


0.224- 


-0.698 


1.37 ±0.35 



Note. - The distances are in units of h^ 1 Mpc. The best fit value r and its 
la confidence level r\ ow — r up are computed from the differential TPCF and the MLE 
method, assuming a fixed value for the slope 7. The values z and z^ are the mean 
redshift of the QSO sample and the median redshift of the observed QSO pairs inside 
20/i _1 Mpc, respectively. The value reported in £(20) is the observed value of the TPCF 
integrated over 20/i _1 Mpc, with its ler confidence level, £i ow ~~ Cup- The bias factor is 
computed assuming the cosmological parameters described in Section 6.1. 



